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Abstract 

The Airy function Ai{z) and its derivative Ai (z) occur in a large 
number of applications in Chemistry and Physics. As a result, there is a 
continuing interest in the properties of these functions. Recently, there has 
been interest in obtaining general expressions for the higher derivatives 
of these functions. In this work, general expressions for the polynomials 
which are contained in these derivatives are given in terms of the partial 
Bell polynomials. 

In view of the continuing occurrence of the Airy functions Ai(z) and Ai (z) 
in a large number of mathematical and physical applications, [1] interest in 
these functions among researchers in Mathematics, Physics and Chemistry has 
provided an impetus for the further exploration of their properties. In this 
paper we investigate the structural properties of a set of polynomials associated 
with the derivatives of the functions Ai{z) and Ai (z). 

The Airy function defined by the equation 

1 f°° 

Ai{z) = - cos(k^ + xt)dt, 

Jo 

with 

Ai{{)) ^ 



32/3r(|)' 



also satisfies the differential equation 

d'^Ai{z) 
dz^ 



zAi{z). 



As a result, direct successive calculations of the higher derivatives of the Airy 
function and of its derivative produces the very simple forms 



d'^Aiiz) 
dz" 

dz^ 



= rn{z)Ai{z)+Qn{z)Ai{z), 

= rn+l{z)Ai{z)+Qn+l{z)Ai{z), 
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where Vn{z) and Q n{z) are polynomials. Higher derivatives of the second 
linearly independent solution of the Airy differential equation, i.e. Bi{z) and 
it derivative Bi'{z) produce exactly the same polynomials. First studied in 
detail by Abramochkin the Abramochkin polynomials can easily be shown 
to satisfy the relations 

V n+2{z) = ZV n{z) +nV n-l{z) , 
Qn+2{z) = zQn{z) +nQn-l{z), 

a z 

Qn+l(.z) = — l-Vniz). 

The three term recurrence equations shown above do not appear to have been 
studied in the past. Nor do the usual methods of solving 3— order difference 
equations appear to be applicable. Generating functions for these polynomials 
follow directly from the recurrence relations and are 

7r[Bi\z)Ai{z + t)-Ai'iz)Bi{z + t)] = V — P„(z), 

n— 

t: [Ai(z) Bi{z + t) - Biiz) Ai{z + t)] = V-rQ„(^)- 

^-^ n'. 

The last two recurrence equations involving the derivatives of Vn{z) and Q „(z) 
are shown in the appendix to be useful in determining the higher derivatives of 
the Airy functions. The first few polynomials are given in the Table I. As seen 
in the table, there are strong regularities in the sequences of these polynomials. 
However, they are also observed to be complicated functions of the index n. As 
a result one would expect that general expressions for these polynomials will 
not be simple. 
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n 


P„(2) 


Q n {z) 


1 





1 


2 


z 





3 


1 


z 


4 


z' 


2 


5 


Az 


z' 


6 


4 + 


6z 


7 


9z^ 


10 + 


8 


28 z + 


12 z^ 


9 


28 + 16 z^ 


52 z + z"* 


10 


100 z^ + z^ 


80 + 20 z^ 


11 


280 z + 25 z"* 


160 z^ + z^ 


12 


280 + 260 z^ + z" 


600 z + 30 z"* 


13 


1380z^ + 36z'' 


880 + 380 z^ + z" 


14 


3640 z + 560 z^ + z*" 


2520 z^ + 4228 z'' 


15 


3640 + 4760z^ + 49z'' 


8680z + 770z* + z' 


16 


22960z^ + 1064z^ + z« 


12320 + 7840z^ + 56z'' 


17 


58240Z + 13160z4 + 64z'^ 


46480z^ + 1400z^ + z« 


18 


58240 + 99120z^ + 1848z'^ + 


151200z + 20160z* + 72z' 



Table I 



Indeed, attempts to obtain general expressions for these polynomials have 
proven to be challenging. In computing the n — derivatives of the Airy functions 
it is useful to note that they are related to the modified Bessel functions of the 
second kind K^{C,) of fractional order v, (sometimes called Basset functions) by 
the relations |3] 



Ai{z) 



A^{z) = -1^/^2/3(0, 
TT V3 



C 



_ 2^3/2 



In the computation of the Vniz) and Q n{z) polynomials it is sufficient to con- 
sider the n — derivative of the function Ai (z). In terms of the Bessel function 
-^2/3(0 these derivatives are given by 



(z) 



d"i^2/3(C) , d"-li^2/3(C) 
z H h n ; 



(1) 



The n — derivative of a function f (Ciz)) which is an implicit function of the 
variable z can be given by Faa di Bruno's formula [4] 
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written in terms of the partial Bell polynomials Bn, k where the polynomials 
B,^. k (of total weight v and degree k in the variables xi, 2:2, ■ • • , Xn^k+i ) are 
given by 

and the sum ranges over all possible values of the integers ki subject to the 
conditions 

ly- k + 1 = I, 

ki + k2 + ■ ■ ■ + ki = k, 
fci + 2 fcs + ^Iki = jy. 

The partial Bell polynomials can easily be obtained using the generating func- 
tion 

I'u-k+l 



.... n\ 

u—k \ /^— 1 



u—k 

from which one gets 

Bn,k{xi,X2,--- ,Xn-k+i) =n\- {[r] — "a*!"^ 



/n-fc+1 ^ 



where the notation ([i"] (X]^=i^ ~/^) ^^^^'^^ ^'^^ coefficient of in 
the sum on the right hand side of the generating function's relation Eq. (3). It 
is useful to note the special cases of Bn.k shown in the Table II. These values 
have been obtain from relations due to Cvijovic |6]. 



k 


Bn,k 








1 




n 




n-1 




n-2 


|(™)xr^[(n-3)xi/xi + |x3] 



Table II 



Alternately, the Bell polynomials can be computed using software packages 
such as Maple (version 15) or Mathematica which contain the partial Bell poly- 
nomials as part of their standard sets of mathematical functions. 
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Using Eq. (2) we obtain the required derivatives as 



(-1) 



ji+i 



7rV3 2« (fC)" 



k=0 



B„,fc(-3!!,-l!!,--- ,(2n-2fc-3)!!) 



-2n ^(30^ 



fe=0 



B„_i,fc(-3!!,-l!!,--- ,(2n-2fc-5)!!)] for n > 1, 



where the k— derivative of C with respect to z as given by 

0=2(-2)-'^(2^-5)!!(|C)^-§^ 

has been used to rescale and simpUfy the required Bell polynomials (the double 
factorial function has values {2k — 5)!! = 1-3-5 • • • (2k — 5) with (—3)!! = — l,and 
(-1)!! = 1). 

The sums above can be combined if the last term in the first sum is written 
out explicitly and use is made of the relation i3„_ „((— 3)!!) = (—1)". Written 
more compactly we have 



d"Ai (z) 



fc=0 



(4) 



where we have defined the term Ai?„^ k as 

AB„, - B„, fe((-3)!!, (-1)!!, • • • , (2n-2fc-3)!!)-2n.B„_i, fe((-3)!!, (-1)!!, • • • , (2n-2fc-5)!!), 

since Xi = (2z — 5)!!. The derivatives of the modified Bessel function in Eq. 
(4) can be obtained from the well known relation 7 for the n — derivative of 
K.iO i.e. 

d^K^iO _ i-ir 



(5) 



Wc note that Bessel functions of order higher and lower than occur in this 
sum. With ^ = 2/3 all of these terms can be reduced to ones which contain 
the functions Ki/^(C,) and i^2/3(C) by means of the relation 



.p=o 



(6) 



-i(C) 



\2p+l) r{ri+p.~p) 



,P=0 



for fi> 0. 



The equation above allows Bessel functions with negative order to be replaced 
with those with positive order. Equation (6) has been obtained by repeated 
application of the recurrence relation [7 for K^{C,) i.e. 
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K^io = if.+2(c) - (f ) {'^ + m^+iio- 

Using Eq. (5) for the derivatives we have 



d'^Ai {z) 



7r\/32" (|C) 



l(n-l) 



n-1 



(-10" E (:)^2/3+2i-„(C) + E(-iO' U E ( ■)^2/3+2i-fe(C) 

(7) 



1=0 



fe=0 



0.1 Calculation of the Bell-free terms of equation (7) 

In order to ease the calculation of the sums in Eq. (7) and without loss of 
generality (cf. the appendix) we consider only even values of n. With n = 2m, for 
m > 1, the resulting calculation oid'^"^Ai {z)/dz'^"^ eventually yields 'P2m+i{z) 
and Q 2m+i (z) . 

Under those eircumstanees the Bell-free term in Eq. (7) which corresponds 
to the even order derivatives of if2/3(C) is given by 

2m m— 1 2m 

E (T)^2/3+2i-2m(C) = (^^')-^2/3(C)+E (T) "^2/3+2 i-2m (0+ E (^r) ^2/3+2 i-2m(C)- 
i=0 i=0 i=m+l 

(8) 

In the first sum of Eq. (8), the orders of the Bessel functions 2/3+2z—2to within 
the range < i < m — 1 are negative and every Bessel function therein can be 
replaced with a term containing the functions if 1/3 (^) and -ftr2/3(C)- Choosing 
v = 2/3 + 2i — 2m, and u + 2n = 2/3 in Eq. (6) requires that ijl = m — i. As 
a result we get 

p=0 

m—i—l 

\^ (m-i-l+p\ r(2/3-m+i+p) p >2p 
+-«-2/3(;j 2_/ \ 2p ) r(2/3-m+i-p) ' 
p=0 

where we have used Eq. (5) and the parity property of the Bessel functions i.e. 
K_^{Q = K^iO for the ii'i/sCC) function. 

In the second sum in Eq. (8) the orders 2/3 + 2« — 2m within the range 
m + 1 < i < 2m are positive. Using the parity property of the Bessel functions 
of order v we get 

■^2/3+2i-2m(C) = ■^-2/3-2i+2m(C)- 

In this case we set v = —2/3 — 2i + 2m, and v + 2fi = —2/3 which gives 
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fi = i — m and using Eq. (6) we raise these indexes to get 

i—m — \ 

V t/-\ — ir ( r\ \ ^ ('i-m— l+p\ r(-2/3-i+m+p) r)^2p 

A_2/3_2i+2ml(,J - 2^ i 2p j r(-2/3-i+m-p) 

p=0 

z— m— 1 

--"■-5/3l.t,J 2^ I 2p+l jr(-2/3-i+m-p)^^ 
p=0 

Employing the relations -?^-2/3(C) = -?^2/3(C) and ii'_5/3(C) = -R'l/sCC) + 
(|)(|) Jr2/3(C)) the latter expression can be written as 



z— m— 1 



-f'^-2/3-2i+2m(C) - --f'^l/3(C) ^ "i^+l) T{--l'l2,-%+m-v) ^C'^^'^ 



p=0 



2 — m — 1 

+i^2/3(C) E 

P=0 L 



(t-m-X+vX r(-2/3-i+m+p) r)^2p |'2^o^2 /^i-m+p^ r(l/3-i+m+p) 

I 2p ^ r(-2/3-i+TO-p)^'' \?,)'^^ 2^ 



p=0 



r(-2/3-i+m-p) 



(10) 

Combining Eqs. (9, 10) and reversing the order of the resulting double sums, 
Eq. (8) can be written in the more compact form 



2m m — 1 

(f E (t)^2/3+2.-2m(c) = -^1/3(0 E (i) 



2rn-2p-l 



i=0 



p=0 



+^2/3(0 



2m in— 1^ 



2rn-2p 



p=0 



or reindcxing the powers of ^ in these sums we have 



do'" E ('r)^2/3+2i-2m(c) = - (f ) i^i/3(c) a^'" E ^^°^("^' ™ - 9) 



i=0 



9=1 



29 



+ i^2/3(C) 



(i|)'"+32"E^2(m,m-g;C)(§) 

9=1 



29 



(11) 



where the coefficient C^^im,^) and the quantity C2(m,p;C) are given by the 
finite sums 



m— 1— p 

r^^^(rr,T,\— \^ r/2m\ /m+p-i\ r(5/3-m+p+») / 2m W»+2p+l\ r(8/3+2p+») 

V'hP) — 2^ [W A 2p+l ^ r(2/3-m-p+j) W+p+l+iA 2p+l ^ r(5/3+i) 
j=0 



and 



m— 1— p 

E 

i=0 



r'.^Cmn-r^- \^ f/2m\ /m-l+p-i\ r(2/3-m+p+») 
<^2im,p,g- 2^ [Uj( 2p ) r(2/3-m-p+i) 



i/ 2m W»+2p+l\ r(8/3+2p+0 r (2p+l) , 2/^2\^11 

"•"Vm+p+l+iA 2p+l / r(5/3+i) 1 (2p+l+i)(5/3+i) ~^ J, \(, ) 
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0.2 Calculation of the Bell-containing terms of equation 
(7) 

In this case, both even and odd order derivatives of the Bessel function i^2/3(C) 
occur. Consequently, when the left hand side of the equation below has been 
rewritten as sums corresponding to the even and odd values of fc, the following 
expression results 

2m— 1 fe m— 1 2 ■ 

(-iC)'=Ai32™,fe5:©i^2/3+2.-fe(C) = E (f ) 'aB2^.2,£(V)^2/3+2,-2,(C) 
fc=0 1=0 j=0 i=0 

m-1 2 1 + 1 

- E (¥) Ai?2„,2, + 1 E CT)K-l/3+2.-2 AO- (12) 

j=0 1=0 

We note that the first sum on the right hand side of Eq. (8) contains a Bell 
term with even degree 2j and the second sum contains a Bell term with odd 
degree 2j + 1. In the relation above it will be convenient to treat the cases with 
J = separately where j has the range < j < rn — 1 . As a result we rewrite 
Eq. (12) as 

2m-l fc 

E (-fC)'=Ai32,„,feE©^2/3+2.-fe(C) = 
fc=0 j=0 

AB2m, 0^^2/3(0 - (f ) AB2m,l [2i^l/3(C) + | (f) ^2/3(0 

+i/(m- 2){E (f)''AS2„,2,E(V)^2/3+2.-2,(C) 

j = l j=0 
m — 1 2j + l 2j + l 

- E (¥) Ai?2™.2, + 1 E f '.+ ')^-l/3+2.-2 ,(C)}, 

j=l 1=0 

where H is the Heaviside [S] step function has been inserted to ensure in the case 
where m = 1 that the double sums don't contribute to the left hand side of the 
equation above. Using the values of AB2m, o = and AB2m, i = — 5(4m — 7)!! 
we get 



E (-iC)'=Ai?2™.feE(')^^2/3+2.-fe(C) = 
fc=0 i=0 

10(4m - 7)!! (f ) [k,^,{0 + (^) ^^2/3(0^ 

, J:T=1 Ai32m,2, E-io (V)^2/3-2,+2.(C) 

+ H [111,-2 ) { 

^ -Er=V(¥) AS2„,2, + lE?io^^e^+^)^-l/3-2,+2.(C) 

The next task being the reduction of the K functions with orders 2/3— 2j' + 2j 
(in the Bell sum with even degree) and — 1/3 — 2j + 2i (in the Bell sum with 
odd degree) to orders 1/3 and 2/3 . 
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0.2.1 Analysis of the Bell sum of equation (12) with even degree 

Wc start with the sum of Bessel functions contained in the Bell sum of even 
degree of Eq. (12) i.e. 

E(V)^2/3+2i-2,(C)=E('/)^2/3+2.-2,(C)+ E f /) ^2/3+2 .-2 , (C) 
i=0 1=0 i=j+l 

(13) 

+ (7)i^2/3(C), for j>l, 

where we have separated the sum of K functions in Eq. (13) into those with 
positive order greater than 2/3, negative order less than 2/3 and order 2/3 re- 
spectively. As seen above in connection with Eq. (6) we have for -K'2/3-2 j+2 j(C) 
in the first sum of Eq. (13) 

i^2/3-2,+2.(c) = -^1/3(0 i^JtrVrm^ (!)''^' 

p=0 

+^2/3(0 r^n ^/S:^ (14) 

p=0 

In the second sum in Eq. (13) we have replaced order 2/3 — 2j + 2i with 
-2/3 + 2j - 2i. We get for K_2/3+2j-2i{0 



j — i — l 



i^-2/3+2,-2.(c)=i^2/3(c) E t%'^nm^^{iy 



p=0 



^-S/SlU 2^ {2p+llT(-2/i-t+j-p)\c) 



p=0 



Using the relation if_5/3(C) = -f^i/3(C) + (2/3)(2/^)i^2/3(C) la^st expression 
becomes 



i^-2/3+2,-2.(C) = -i^V3(C) E' C^p^) :^^^) (!)''^' 

p=0 

, (r\ r/i-i-l+p\ r(-2/3-i+j+p) 2 (2\ (i-j+p\ m/3-z+j+p) 1 /'2\ ^ 

+ -«-2/3ig 2^ [( 2p jr(-2/3-i+j-p)-3 I^CJ l2p+l jr(-2/3-i+j-p)i [cj ' 
p=0 

(15) 

Using Eq. (14) and Eq. (15) we get after reversing the order of summation in 
Eq. (13) a relation in complete analogy to the result found in connection with 
the Bell free term in Eq. (11) i.e. 
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2 J j \ j 1 

1=0 p=0 p=0 

(16) 

As a result, the Bell term of even degree in Eq. (12) can be written as 

rn — 1 ■ 2 7 



E (f) 'AB2r„,2,X(V)^2/3-2,+2i(C) 



m— 1 / m— 1 



2? 



-(f )^i/3(c) E E 3'^' ^^^^.^i • - 9) (!) 

{m — 1 / m — 1 \ 2 m— 1 2 ■ 

^ K]32^Ai?2„,2, •C2(j, j-g;C) (f) ' + E (i ) ^^2^.2, (^^ ' 

(17) 

0.2.2 Analysis of the Bell sum of equation (12) with odd degree 

We now consider the sum of Bessel functions contained in the Bell sum with odd 

degree of Eq. (12). As before, we have separated the sum into Bessel functions 
with order 1/3, and those with positive order greater than and less than 1/3 i.e. 

2j+l j-l 2j+l 

Y: f^+^)^-l/3-2,+2.(C)=E(''r^)^-V3-2i+2.(C)+ E f^+^)^l/3+2,-2.(C) 
i=0 1=0 i=j+l 

+ C^+')Kys{0, for j>0. 
Using Eq. (6) we obtain for j > 0, 

, (n - K , \^ rp-l-»+P\ r(-l/3-j+»+p) 1(2^ (3-i+P\ r(2/3-j+i+p) ^ ( 2\ 

A_i/3_2j+2i(.g - Ai/3^g 2^ [(^ 2p J r(-l/3-j+i-p) ~3 ) \ 2p+l ) T{-l/3-j+i-p) \ \Z ) 



-(1)^2/3(0 E T{-y3-t+t-p) (f) ' 



p=0 



•'-I . . 2p 



^^1/3+2,-2.(0= ^1/3(0 E t%'^rT%^^{D 

p=0 

2— 1 „ 

-[^)^2/3[Q 2^ l2p+i jr(i/3-i+j-p) [c) ■ 



p=0 
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Using these results we get for the Bell sum of odd degree 

2j + l 



=0 

+E 

p=0 



) ST I'^J+l'l (j-i+P\ r(4/3+j-z+p) / L 1 f 2 A 

A i A2p+i;r(l/3+j-i-p) I 0-i+p)(l/3+j-i-p) 3 ^^Cy* 

(!) ^^/3(C)[E (e G.^.UCtr.r)^^¥^| 



r(5/3+»+2p) 

r(5/3+i) 



p=0 



E(!) E m( 



p=0 



j-i+p\ r(4/3+j-»+p) 
2p+i J r(i/3+j-i-p) 



], for j>0. (18) 



i=0 



In the equation above we have interchanged the order of summation in the 
double sums and have replaced the ratios of Gamma functions with nega- 
tive arguments with ones with positive arguments using the relation T{—z) = 
— TT csc(7rz)/r(z + 1). The expression in Eq. (18) can be simplified by writing 
out the p = terms which occur in the double sums whose ranges are < p < j 
and then re-indexing those sums. We get as a final expression 



2j+l 



E r^+')^-l/3-2 ,+2.(0 = ^1/3(0 



j=0 



2^^ + (T)+E^3(i,p;C)( 

p=0 



2p 



(f) Ky.iO (i + 1)2'^ + \{j + l)C^I')+Y,C,{j,p;0 (I)'' 

p=0 



, (19) 



where 



csu,p;o= E imaryrmi^ 



i=0 



(l+2p) 



0+p-j)(l/3+j-i-p) ^ 3 



+ 



,( 2j+l W»+2p+2\ r(n73+2p+i) /2A^ 1 f • ^ 
+ [i+j+p+2){ 2p+2 ) r(5/3+i) [c) >' tor J > U, 



and 



i-i-p 



r-rin-^^- \ " r/^2j + lWj+p-i\ r(4/';i+j-i+p) / 2j + l Wi+2p+3\ r(J_l73+2£+4 
L'4U)/^!'5J— 2^ 11 i J\2p+1 Jr{l/3+j-t-p)^\i+j+p+2j\ 2p+3 J r(2/3+i) 



i=0 



i(f)'}, for j>0. 



As a result, the Bell term with odd degree in Eq. (12) can be written as 



m— 1 

E 



2i+l 



2j+l 



AB2m,2j+l Cf')K-l/S+2i-2j{0 



i=0 
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(!) ^i/3(c)[ E {2^^ + rr)}3^^-+^AB,„,.,+i ( 

m-1 / m-1 \ / s 2fe 

+ E E 3'^"+^AB2rr>,2,+lC30-,J - fc; C) (§) ] 
fc=l \j=fc / 

ni—1 



2j 



m— 1 / m— 1 



+ E E 3^^+^AB2™,2,+iC74(j,i - fc; C) (§) ] 



(20) 



k—l \ J — A; 



Where we note that in the case where m = 1 the sums in Eq. (20) are empty. 
Gathering the results in Eqs. (12, 17, 20) we have for 



1 



7r\/322'" (|C) 



3 .s|(2m-l) 



where 



a{m; Q = i^i/3(C)[30 (4m - 7)!!(|) - 3^- E ^(™' ^ " ^ " 



j=0 



m— 1 



- Him - 2) E [2^^' + CT)] ^B2m,2j+iC4r'^' 



if(m-2) E E {3'^AB2„,2,Cl°'(i,i-g)(i)2''-i+32^+iAB2^,2,-+iC3(i,i-g;C)}] 



g=l j=g 

m 

+i^2/3(C)[10(4m - 7)!! + {'Z') + E '^^(m, m - j; C)(|)'^' 

m— 1 

+ if(m-2) E [(^■)Afi2„.,2, - {(1 + Ij)2'^ + (j + 1)(%+^)} AiJ2„,2,+i] (f ) 

m—1 m—l 

+ H{m-2) E (I)'' E 3'^[Ai32„,2, C2(j,j-(?;C)-3Ai?2„,2,+iC4(j, j-g;C)]]. 

9=1 j=9 

(21) 

In the equation below we have combined terms appearing in Eq. (21) noting 
that AB2m,i = — 5(4m — 7)!! and having rewritten the sums in preparation for 
the extraction of the sought-after polynomials from that expression. We have 



2j 



(7(m;C) = - 



m — l 

KiMOi E {3'" )(m,m-l-j)+[22i + (2i+i)]32i+iAi32™,2,+i}(i)'^^ 
3=0 
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+H{m-2) J2 E ^''{c[°\j,J-q)AB2„,,2j+i iC/2f ■Cs{j,j-q;0 AB2rn,2j+i}] 

q=l j=q 

m 

+i^2/3(C)[10(4m - 7)!! + {^Z) + 3'" E m - j; C)(|)'^' 

m— 1 

+ E [(i')^^2„.,2, - {(1 + W + iJ + )} Ai?,„,,,+i] (f )^^- 

rn — 1 m — 1 

+ iJ(m-2) Ed)'' E 3''{^2(j,i-g;0 AS2„x,2j-3C4(j,j-g;C) AB2™,2,+i}]. 

9=1 j=9 

(22) 

We have 

-7, ^ io-i/3(™;C) ^ . . cr2/3(n^;C) 



22m (3 ^ / 2)''(™-^)/^ ' ™+ V y (3 C / 2)^"^^ ' 

where 



m— 1 



<^i/3(m;C) = E {3'"Cf'Km-l-i) + [22^ + f^+i)]32^+iAB2™,2,+i}(|)'^' 

m— 1 m— 1 

+if(m-2) ^ 32^"{cj"'0-, j-g) AB2m,2i + 3 (C/2)'-C3(j, j-g; C) Ai?2„,2,+i}, 

9=1 j=9 

(23) 

and 

rn 

^2/3(m; = 10(4m - 7)!! + (t') (f )'™ + 3^™ E ^^(m, m-j; C)(i)^^" 

j=i 

m— 1 



+ if(m-2) Y: [(^■)AB2„,2, - {(1 + ii)2^^- + (i + 1)(%+^)} AB2™,2,+i] (f 



j=i 

m— 1 m— 1 

+ if(m-2) E (i)'" E 3'''{C2(i, j-g; C) AB2m,2j-3C4(i, j-q; Q AB2m,2j+i}. 
9=1 j=q 

(24) 

In the expressions for ■P2m+i(-z) and Q2m+i{z) given above, cancellations 
of terms with negative powers of ( must occur in ai/^/K / 2)'^'""^ ^■'^'^ and 

f2/3/(C/2)^™^'^ leaving terms with only zero or positive exponents as is re- 
quired in the sought-after polynomials. To this end, we rewrite the terms in 
c7i/3 and (T2/3 gathering together ( terms of the same order. In this regard it 
is helpful to write the coefficients C2, C3, C4, as 

C2{m,p;C) = (m,p) + {2/CfC^^\m,p), 
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C3(m,p;C) = (m,p) + (2/C)'cf 
C4{m,p;0 = Cf (m,p) + (2/C)'cf (m,p). 
In the case of to > 4, we have 



(C/2)4(™-i)/3 



m-2 



^ [A{m,q)+D^^\m,q)'] (C/2)2[9-2(™-i)/3l + ^ i^W (m, g+l)(C/2)2[«- 

9=2 



2(m-l)/3] 



9=2 



where 



A{m,j) = 32™c(°)(to,to- 1 - J-) + [22J- + (2^+i)]32j"+iAB2,„,2j+i, 

m — 1 

i?(0)(TO, q)=Y. 3''{Cf J - g) AB2m.2j + 3Cf J - q) Ai?2™,2, + l}, 



J=9 



i?(2)(TO,,) = ^32^+icf (j,j-g)Ai?: 



2m,2j + l- 



Since the exponent q — 2(to — 1)/3 > the sums above may be replaced by 

m — 1 m — 2 



o"i/3("^;C) 

(C/2)4(™-l)/3 



^ [A(TO,g)+i?(2)(„^^^)j (^/2)2[«-2('"-i)/3l+ ^ i?(")(TO,g+l)(C/2)2['- 

q=M(m) 

(25) 



2(m-l)/3] 



r;=M(' 



where 



M(to) 



2 (m-1) 
3 



and [a:] is the ceihng function whose value is equal to the smallest integer 
which is greater that or equal to a;. A plot of the function M(to) is shown in 
Figure [TJ 

Reindexing the sums in Eq. (25) we have 



gi/3("i;C) 

(C/2)4('"-l)/3 



M(m) 



^ \A{m,q + U{ra))+D''^\m,q+U{m)) (C/2)2['? 



+ M(r?x)-2(m-l)/3] 



9=0 



M{m) -1 

- ^ i?(0)(TO,g+M(TO) + l)(C/2)2[9+^(™)-2(m-l)/3]^ 
9=0 



where 



7\/((to) = to — 1 — M(to). 
A plot of M (to) is seen to be staircase function of m as shown in Figure [21 



14 



0.0 



2.0 



4.0 



6.0 

m 



8.0 



10.0 



12.0 



Figure 1: 



Using the result in Eq. (23) we get (with ( = ^z^/"^) 



2m+l 



, e(m) 



2217132 M{r, 



-j-{ A(m,m — 1) + D'-^\m,m-l) \ (zV9) 



M{m) -1 

^ Cv3(m,g) (zV9)'}, 

9=0 



where 



Ci/3(m, g) = A{m,, M(m) + q) + D'^°'>{m, M(m) + 1 + q) + D^^\m, M(m) + g), 
and where 

e(m) = 3 M(r7i) - 2 (to - 1). 

It should be noted that the exponent e(m) is a periodic function of m i.e. 
e(m) = e{m + 4) as seen in Figure |3l 

0.3 General expressions for V2m+i{z) and Q2m+i{z) 

The expression in Eq. (26) can be further simplified by evaluating the sums 

which are contained in the coefficient of in conjunction with the 
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Figure 2: 

relation Ai?2m,2m-i = 2to^ + 3m. We get 

M(m) -1 
3=0 

where the coefficients Ci/3(m, g) are given by 

o2M{m) + l r , , 

Ci/3{m,q)= [c[°\m,M{m)-q) + m{l + 2m/3)C^°\m-l,M{m)-q) 



with 



rr,2q-2M(m)-2 , o-2m /'2M(rri)+l+2q\ -i ad , q 

12" ^ ^ +2 ( M(„)+, j|A-»2m,2M(m) + l+2g +2 



M{m)—l—q 

-2m, ^ Ci/3 (m,g,i) 

3=0 

(27) 



,2M(m) + l+2g+2j 

+3 C[°^ {M{m) + l + q + j,j) AB^m ,2M(m)+2+2g+2j 
+9 Cf ^ (M(m) + l + q + j, j) AB^m ,2M(m)+3+2q+2j- 

and where we note that the sum in Eq. (26) contributes to V2m+i{z) only if 
m > 4. 
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Figure 3: 

In a similar way, the polynomials Q 2m+i {z) can after some simplification 
be written as 

A4(m+l)-2 

Q2m+i iz) =z^+ (™) (3m + 1) ^"-3 + ^^(-+1) . C2/3(m, q) (^3/9) ' , 

9=0 



where the coefficients C2/-T,(m, q) are given by 



C2/3(m,g) = " [cf\mM{m+l) - q) + C^^\m,A1(m + 1) -l-q)] 



(2), 



629 



+ 22-M('"+l) [^^2m,2M(m+l)+2g - {2 + 3M(m + l)+3g'} i^B^mfiM {m+l)+l+1q] 



329+2 



M(m)-2-q 



i=o 



with 



'C2/3{m,q,j) = 

= + M(m + 1) + 1 + g, j + 1) + Cf\j + M(m + 1) + 1 + j)\ ASa^, 2M(m+i)+2+29+2j 

-3 |cf (j + M(m + 1) + 1 + g, j + 1) + Cf\j + M(m + 1) + 1 + g, j)| ABa™, 2M(m+l)+3+2g+2j) 
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and where we note that the sum over C2/3(m, g) is empty for m < 5 and the 

relation 

3 [AB2m, 2m -2 - (3m - 1) ABsm, 2m -1 + (27to - 7)] = C^) (3m + 1), 

has been used to simplify the coefficient of z"^~^ in the expression for Q2m+iiz)- 

It does not appear to be possible to further simplify either of the Ci/3(m, q), C2/3(m, q) 
coefficients. In practical terms, the expressions for V2m+i{z) and Q2m+i{z) are 
seen to exhibit a structure which depends on the discontinuous functions M (m) , 
Ai{m) and e(m) and as a result these functions are in practical terms responsible 
for the structural nature of those polynomial sequences. Ultimately however, 
the underlying structure and properties of these polynomials rests upon the 
nature of the Airy function itself. 

Appendix 1 



Having obtained general expressions for the polynomials of odd order i.e. 
'P2m+i{z) and Q2m+i (^) , the even order polynomials can be obtained from 
the recurrence relations given above. We have 

V2m{z) = [P2m+3{Z) ' zV2m+M] I i'^rU + 1), 
Q 2m {z) = [Q 2m+3 {z)-zQ 2m+l {z)] /(2m + 1). 

Alternately the even order polynomials can also be obtained using the differen- 
tial difference equations also shown above i.e. 

d'P2m-l{z) 



'P2m{z) = zQ2m-l (z) + 
Q2m {z) = V2m-l{z) + 



dz 

dQ 2m-l {z) 



dz 

Appendix 2 

Expressions for the coefficients C^'^\m,p) and Cj^^'(m,p) have been given 
below. 



m—l—p 



r^(°Vrr) - \^ [ (2m\ (m+p-i\ r{5/3-m+p+i) I 2m \ (i+2p+V\ r(8/3+2p+») 
— 2^ \\ i ) \ 2p+l y r(2/3-m-p+i) W+p+l+iA 2p+l ) r(5/3+i) 

i=0 

m—l—p 

r-C^) /'rr, T,^ - ^/'2m^ /-m-l+p-i\ r(2/3-m+p+i) , ( 2m \ (i+2p\ r(8/3+2p+j) i 

<^2 y^^P) - 2^ [\ i )\ 2p jrW3^^^^^^=P+iy + (m+p+l+J( 2p jTW3+i]~J 
i=0 

m—l—p 

M'^)(rr, r,^ - 1 \^ ( ^rn \ /»+2p+l\ r(8/3+2p+i) 
^2 V'hP)—^, Krn+p+l+i)\ 2p+l I r(5/3+i) 
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Cf\mp)^ ^ (^^m+lj (^n+p-l-i\ r(4/3+m+p-») 



2p / r(4/3+m-p-j) 
i=0 



m— 1— p 



r^'^)(mr,\— \^ n arn+l\(m+p-i\ T(4./?,+m+p-i) , ( 2m+l Wi+2p+2\ r(n/3+2p+i) t 
•-^S V'hP] — 2^ i A 2p+l yr(l/3+rn-p-j)+Vj+m+p+2A 2p+2 / r(5/3+i) J 



m — 1 - 



1^0 

711— 1—p 

r'^^^(m'n\- \^ ( 'I A+2P+3') r(ll/3+2p+i) 

O4 \in,p)~ Vj+m+p+2A 2p+3 / r(2/3+j) 

The finite sums which make up these coefficients can also be represented by 
combinations of the hypergeometric functions 3^2 and 4F3 although there ap- 
pears to be little or no advantage gained in expressing them in that way. In 
less general terms each of these sums can be written as weighted polynomials in 
m. For example cf^ (ra^j)) = 2^™(^^^_j^)-p(P-'(TO) where the symbols p^'^\m) 

or p'*^")(m) indicates a polynomial in m of order uj. In the case of the other 
coefficients one finds 

C«(m,p) = 22'"p(2p+i)(m) + p'(2p+i)(^).(2™), 



V m+1 / ' 



Cf^ {m, p) = 22'"p(2p) (m) + p' (m) 
C^\m,p) = 22'"p(^''+i)(m) + p'(2^+i)(m).f:V+2) 
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